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Abstract 

A fully self-consistent calculation of the bosonic dynamics of the Hubbard model is developed within the Composite 
Operator Method. From one side we consider a basic set of fermionic composite operators (Hubbard fields) and 
calculate the retarded propagators. On the other side we consider a basic set of bosonic composite operators (charge, 
spin and pair) and calculate the causal propagators. The equations for the Green's functions (GF) (retarded and 
causal), studied in the polar approximation, are coupled and depend on a set of parameters not determined by the 
dynamics. First, the pair sector is self-consistently solved together with the fermionic one and the zero-frequency 
constants (ZFC) are calculated not assuming the ergodic value, but fixing the representation of the GF in such a 
way to maintain the constrains required by the algebra of the composite fields. Then, the scheme to compute the 
charge and spin sectors, ZFCs included, is given in terms of the fermionic and pair correlators. 
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Recently, the Green's function method for compos- 
ite operators has been revisited [1] . In particular, it has 
been shown that the formulation generates an internal 
self-consistency which cannot be solved uniquely by the 
dynamics, but ingredients related to the microscopic 
nature of the local operator algebra and to the macro- 
scopic nature of the external boundary conditions must 
be provided. This is not surprising. The properties of 
composite operators are not known at priori, they have 
a microscopic nature but manifest at level of observa- 
tion; as a consequence they are self-consistently deter- 
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mined by the dynamics of the system, by the algebra 
and by the boundary conditions. 

Roughly, the properties of electronic systems can be 
classified in two large classes: single particle proper- 
ties, described in terms of fermionic propagators, and 
response functions, described in terms of bosonic prop- 
agators. These two sectors, fermionic and bosonic, are 
not independent, and a fully self-consistent solution 
requires that both sectors are simultaneously solved. 
All the new theoretical schemes, developed in the last 
years, show the importance of the spin and charge cor- 
relations in order to describe the single particle prop- 
erties of highly correlated electron systems. 

In order to illustrate these ideas, we consider the 
Hubbard model, described by the Hamiltonian 

H = ^(tij-At5ij)ct(i,t)c(j,t) + (/^nT(i)n|(i) (1) 

We use the standard notation: c{i), (i) are annihila- 
tion and creation operators of electrons in the spinor 
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notation; i stays for the lattice vector and i — 
fj, is the chemical potential; denotes the transfer in- 
tegral; U is the screened Coulomb potential; na{i) = 
cJ-(i)co-(j) is the charge density of electrons at the site 
i with spin a. For a cubic lattice and by considering 
only nearest neighbor sites the hopping matrix has the 
form tij = —2dtaij, where d is the dimension and aij 
is the projection operator 



k 71=1 



We choose as fermionic basis 




(2) 



(3) 



where ^(i) = [1 — n(i)]c(i) and "qii) = n{i)c{i) are 
the Hubbard operators, and n{i) — X^^7iCT(i)- In the 
two-pole approximation [2] the retarded GF G{i,j) = 
{R[ip{i)'4>^{j)]) satisfies the equation 

[a;-e(k)]G(fc,a>)=/(k) (4) 

where 7(k) = F.T.{{iP{i,t),ilj^ {j,t)}) and e(k) = 



F.T.{{ 



8jp(i,t) 

at 



iP^(j,t)})I-\k.); the symbol F.T. de- 



notes the Fourier transform. In the paramagnetic phase 
the energy matrix £(k) depends on the following sot of 
internal parameters: H, A = (C"(«)C^('0)^{^"(*)^^(*))> 
P = (n«(i)n^(i))/4- {[cT(i)c|(i)]Vt(i)ct(z)), which 
must be self-consistently determined. Given an opera- 
tor ^(i), we are using the notation ^"(i) = OijCUi*)- 
The operator n^{i) = c^(i)crpc(i) [cr^ = (l,cr), a are 
the Pauli matrices] is the charge {fj, = 0) and spin 
(/i = 1, 2, 3) density operator. The local algebra satis- 
fied by the fermionic field (3) imposes the constraint 
(^(i)?7^(i)) = 0: this equation allows us to solve self- 
consistently the fermionic sector. However, it is worth 
noticing that the presence of the parameter p directly 
relates the fermionic sector to bosonic sectors. 
We consider then the composite bosonic field 



(5) 



where p^(i) = {i)aij,c°'{i) — c™^(i)cr;^c(i). In the two- 
pole approximation [3] the causal GF G^'^\i,j) = 
{T[N^^'\i)N'^''^\j)]) satisfies the equation 

[u;-£'^'(k)]G<''>(k,w) = (k) (6) 

where /('''(k) = F.T.{[N^>'\\,t), N'^>''>\S,t)]) and 
^ F.r.([i ^^'g/'-'\ iV<^)t(j^^)])[j(M)(;^)]-i^ In 

the one-dimensional case (we consider the ID case 
just for the sake of simplicity) the energy matrix 
e''^''(k) depends on the following set of internal pa- 
rameters: (i) fermionic parameters: C™ = {c" {i)c\i)) , 




Fig. 1. Pair energy spectra a?^^^(fc) as a function of the mo- 
mentum k for Coulomb repulsion U = i, temperature T = 0.01 
and filling n = 0.8, 0.9 and 1. 
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Fig. 2. Pair zero-frequency constant P^^^ (0) as a function of 
the filling n for temperature T = 0.01 and Coulomb repulsion 
U = 1,2 and 4. 

C = {c^{i)c\i)), = (c^(i)ct(i)), E = mrH:^)), 
E^ = {'?''(*) c^(*))i where rjij, \ij are the projection 
operators on the second and third nearest neigh- 
bors, respectively; (ii) bosonic parameters: a^, bp 
and c^, whose explicit expressions, although for the 
2D case, are reported in Ref. [3]; (iii) zero-frequency 
functions (ZFF) r(f''(i,j) (sec Ref. [1]). Due to the 
hydrodynamic constraints, two bosonic parameters 
should be determined as = a^i + n — 2{E — E^) 
and = — n + 2{E — E^). The parameter in- 
stead can be determined by means of the local algebra 
constrain {n^{i)n^{i)) = (n) + 2(25^.0 — 1)-D, where 
D = (n)/2 — E is the double occupancy. The ZFF are 
left undetermined. 

According to this, we need another composite 
bosonic field 



B{i) 



where h{i) = C|(i)c|(i) and f{i) = cf(i)c|(i) -|- 
C|(i)cJ(z). In the two-pole approximation the causal 
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GF G^^\i,j) = {T[B{i)B\j)]) satisfies the equation 
[u - e*''' (fc)]G<''' (fc, = (fc) 

where /(^'(fc) = t), fity^ j)]) ande(-s)(fc) = 
F.r.([i^^i^, The energy matrix 

e*-^^ (fc) depends on the following set of internal param- 
eters: (i) fermionic parameters: C", C, C^, i?, S''; 
(ii) bosonic parameters: u, v and w, whose lengthy ex- 
pressions are not reported here for the sake of brevity; 
(in) ZFF r<^'(i,j). The condition that the pair en- 
ergy spectra Lul^\k) are finite and the local algebra 
constrain {h{i)h^ (i)) = 1 — (n) + D completely de- 
termine the parameters u, v and w. The algebra con- 
straint {h{i)f'' (i)) = 2{^°'{i)c^ (i)) will be used to com- 
pute the ZFC r'f (0). Then, the pair sector can be 
immediately solved once the solution for the fermionic 
sector has been found. The results of this scheme are 
shown in Figs. 1 and 2, where the pair energy spectra 
and r5^^'(0) are shown, respectively. It is worth notic- 
ing that, in the present scheme, the pair dynamics is 
ergodic only at half-filling. 

Now, once we have solved the fermionic and pair 
sectors we can come back to the charge-spin one. We 
need to compute the six ZFC: r^°' (0), Pji' (a), r[\^ (0), 

Tn (a), r^2'(0), r^2'(0)- They can be fixed by means 
of as many algebra constrains coming from the ex- 
pressions of the following correlators: {n^(i)n^(i)), 
(Pm(»)Pm(»)>. (h'^iiWii)) and The explicit 

expressions of the constraints are quite lengthy and 
will be given elsewhere. 

In conclusion, we have reported a fully self-consistent 
scheme of calculations for both the fermionic and the 
three (spin, charge and pair) bosonic sectors of the 
Hubbard model. It is worth noticing that, within this 
scheme, the ZFC of the spin and charge sectors, which 
could assume, at least in principle, not ergodic values 
as the pair one does, can be self-consistently computed 
and give invaluable information regarding the dynam- 
ics of the corresponding operators. 
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